Symmetry operators for Dirac's equation on two-dimensional spin manifolds * 



Lorenzo Fatibene^ 
Dipartimento di Matematica, Universita di Torino, Italy 
IN FN - Iniz. Spec. Nal2 

Raymond G. McLenaghan 
Department of Applied Mathematics, University of Waterloo, Waterloo, Ontario, N2L 3G1, Canada 



oo 

o 

O 

(N 

O 

Q 



Giovanni Rastellp 
Dipartimento di Matematica, Universita di Torino, Italy 

Shane N. Smith^ 

It is shown that the second order symmetry operators for the Dirac equation on a general two- 
dimensional spin manifold may be expressed in terms of Killing vectors and valence two Killing 
tensors. The role of these operators in the theory of separation of variables for the Dirac equation 
is studied. 

PACS numbers: 04.20.Gz, 02.40.Vh, 04.20.q 



Oh- 



> 
On 

(N 
(N 

(N 
OO 

o 



I. INTRODUCTION 



Exact solutions to Dirac's relativistic wave equation by 
means of the method of separation of variables have been 
studied since the equation was postulated in 1928. In- 
deed, the solution for the hydrogen atom may be obtained 
by this method. While there is a well developed theory of 
separation of variables for the Hamilton-Jacobi equation, 
and the Shrodinger equation based on the existence of va- 
lence two characteristic Killing tensors which define re- 
spectively quadratic first integrals and second order sym- 
metry operators for these equations (see [1, 10, 13, 18]) 
an analogous theory for the Dirac equation is still in its 
early stages. The complications arise from the fact that 
one is dealing with a system of first order partial dif- 
ferential equations whose derivation from the invariant 
Dirac equation depend not only on the choice of coor- 
dinate system but also on the choice of an orthonormal 
moving frame and representation for the Dirac matri- 
ces with respect to which the components of the un- 
known spinor are defined. Further complications arise 
if the background space-time is assumed to be curved. 
Much of the progress in the theory has been stimulated 
by developments in Einstein's general theory of relativity 
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where one studies first quantized relativistic electrons on 
curved background space-times of physical interest such 
as the Schwarzschild and Kerr black hole space-times. 
This work required the preliminary development of a the- 
ory of spinors on general pseudo-Riemmanian manifolds 
(see [7], [5], [8], [9]). The solution of the Dirac equation in 
the Reissner-Nordstrom solution was apparently first ob- 
tained by Brill and Wheeler in 1957 [2] who separated the 
equations for the spinor components in standard orthog- 
onal Schwarzschild coordinates with respect to a mov- 
ing frame adapted to the coordinate curves. A compara- 
ble separable solution in the Kerr solution was found by 
Chandrasekhar in 1976 [6] by use of an ingenious sepa- 
ration ansatz involving Boyer-Lindquist coordinates and 
the Kinnersley tetrad. The separability property was 
characterized invariantly by Carter and McLenaghan [4] 
in terms of a first order differential operator constructed 
from the valence two Yano-Killing tensor that exists in 
the Kerr solution, that commutes with the Dirac opera- 
tor and that admits the separable solutions as cigenfunc- 
tions with the separation constant as eigenvalue. Study 
of this example led Miller [19] to propose the theory of a 
factorizable system for first order systems of Dirac type 
in the context of which the separabilty property may be 
characterized by the existence of a certain system of com- 
muting first order symmetry operators. While this theory 
includes the Dirac equation on the Kerr solution and its 
generalizations [12] it is apparently not complete since 
as is shown by Fels and Kamran [12] there exist systems 
of the Dirac type whose separablity is characterized by 
second order symmetry operators. The work begun by 
these authors has been continued by Smith [20], Fati- 
bene, McLenaghan and Smith [11] and McLenaghan and 
Rastelli [15] who studied the problem in the simplest 
possible setting namely on two-dimensional Ricmanian 
spin manifolds. The motivation for working in the low- 
est permitted dimension is that it is possible to examine 
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in detail the different possible scenarios that arise from 
the separation ansatz and the imposition of the separa- 
tion paradigm that the separation be characterized by a 
symmetry operator admitting the separable solutions as 
cigenfunctions. The insight obtained from this approach 
may help suggest an approach to take for the construc- 
tion of a general separability theory for Dirac type equa- 
tions. Indeed in [15] systems of two first order linear par- 
tial equations of Dirac type which admit multiplicative 
separation of variables in some arbitrary coordinate sys- 
tem and whose separation constants are associated with 
commuting diffrcntial operators are exhaustively charac- 
terized. The requirement that the original system arises 
from the Dirac equation on some two-dimensional Ric- 
mannian spin manifold allows the local characterization 
of the orthonormal frames and metrics admitting separa- 
tion of variables for the equation and the determination 
of the symmetry operators associated to the separation. 

The purpose of the present paper is to take this re- 
search in a different but closely related direction. Fol- 
lowing earlier work of McLenaghan and Spindel [17] and 
Kamran and McLenaghan [14] where the first order sym- 
metry operators of the Dirac equation where computed 
on four-dimensional Lorenzian spin manifolds and McLe- 
naghan, Smith and Walker [16] where the second-order 
operators were determined in terms of a two-component 
spinor formalism, we obtain the most general second- 
order linear differential operator which commutes with 
the Dirac operator on a general two-dimensional Rieman- 
nian spin manifold and show that it is characterized in 
terms of Killing vectors and valence two Killing tensors 
defined on the background manifold. The derivation is 
manifestly covariant; we work in a general orthonormal 
frame and do not chose a particular set of Dirac matrices. 
The paper is organized as follows: In Section 2 we derive 
the conditions that must be satisfied by the coefficients 
of a second order linear operator which commutes with 
the Dirac operator. In Section 3 we expand each zero- 
order coefficient operator in terms of a suitable basis of 
the Clifford algebra C(2) and solve sequentially the three 
defining conditions obtained in Section 2. By this pro- 
cedure we show that the general second order symmetry 
operator is characterized by a valence two Killing tensor 
field, two Killing vector fields and two scalar fields de- 
fined on the background spin-manifold. These fields are 
related by a certain differential condition the intcgrabil- 
ity condition of which is solved on the Liouville surface 
which is the most general surface which admits a general 
valence two Killing tensor. In Section 4 we give theorems 
characterizing non-trivial first order symmetry operators 
in terms of Killing vectors and non-trivial second order 
symmetry operators in terms of characteristic valence two 
Killing tensors. We also give an explicit example of an 
non-trivial, irreducible Killing tensor which is defined on 
the asymmetrical ellipsoid surface. In Section 5 we estab- 
lish a link between the second-order symmetry operators 
studied in the previous sections and the naive separation 
of variables considered in [15] by exploiting the fact that 



the existence of a non-trivial second order symmetry op- 
erator implies in general that the Liouville metric has at 
least one ignorable coordinate. Section 6 contains the 
conclusion. The notation and conventions of this paper 
are consistent with [15]. 



II. TWO-DIMENSIONAL SPIN MANIFOLDS 

We fix Euclidean signature by setting r] = (2, 0) and 
denote by rj a (, = S a b the corresponding canonical sym- 
metric tensor. A representation of the Clifford algebra 
C(2) is induced by a set of Dirac matrices j a such that 

7a 76 + Ibla = 2»? o6 I (1) 

We stress that we shall never fix a particular set of 
Dirac matrices, but we shall use only algebraic conse- 
quences of (1). The even Clifford algebra is spanned by 
I, 7 a , and 7 := 7172 and the most general element of the 
corresponding Spin(2) group is 

S = al + 67 with a 2 + b 2 = 1 (2) 

One can define a covering map I : Spin(2) — > SO(2) by 
setting 

We stress that until now we are at a purely alge- 
braic level. Let M be a connected, paracompact, 2 di- 
mensional spin manifold. Let P — * M be a suitable 
Spin(2)-principal bundle, such that it allows global maps 
A : P — > L(M) of the spin bundle into the general frame 
bundle. The local expression of such maps is given by 
spin frames e^; see [7]. A spin frame induces a metric 
9nv — r\ab e-1 and a spin connection 

T$ = e%{T a Pll e hf) + d l >e ha ) (4) 

where r Q ^ denotes the Levi-Civita connection of the 
induced metric t^„. We note that such a connection 
satisfies V^e^ = 0. We also remark that Latin indices 
are raised and lowered by the inner product r\ a b while 
Greek indices are raised and lowered by the induced met- 
ric g^ v . For subsequent use we introduce the frame co- 
variant derivative V a := e^V M . The Dirac equations 
then has the form 

Dip = ij a V a ip -mip = 0, (5) 

where the covariant derivative of the spinor tp is defined 

as 

V^ = d^ + iea r; b 7 ^ (6) 
and e a b denotes the permutation symbol. 
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A spin transformation is an automorphism of the spin 
bundle P, i.e. locally 

(*' = /(») 

\g' = S(x)-g {) 

They form a group denoted by Aut(P) which acts on 
spinors, frame and spin connection as 

y/ = s • v 

e'S - ^e^(5) (8) 

r'f - (r^(5) + d^(S) v db ) 

leaving the Dirac equation (5) invariant. Here is the 
I 



Jacobian matrix of the spacetime transfomation x' = 
f(x) onto which the spin transformation projects and 
denotes the anti-Jacobian. 
A second order symmetry operator for the Dirac equa- 
tion is an operator of the form 

IK = E ab V afc + F a V a + GI (9) 

which commutes with the Dirac operator D. Here V ab = 
\ (V a Vb + V;,V a ) denotes the symmetrized second co- 
variant derivative (expressed in the frame). The coeffi- 
cients E ab , F a , G are matrix zero-order operators. 

The symmetry equation [K, D] = 0, can be easily ex- 
panded to the equivalent equations for the coefficients 



'E^V _ 7 ( c E a6 ' = 

f(ajb) — ^( b ][r a ) — <-y c V c E ah 

' G7 a - 7 a G = 7 c V c F a - f (E a V + 7 c E ab ) e bcl + f (E b <V + 2 7 C E M ) e a d e bc 
j a \7 a G = f (F a 7 fc + 7 b F a ) 7 e ab + ^ (2E afc 7 c + 7 c E ab ) 7 e ac V fc P 



Here R denotes the Ricci scalar of the induced metric. 
We observe that, in 2-dimensions, it encodes all the cur- 
vature information of the surface. In particular the frame 
components of the Riemann tensor may be written as 



R a 



bed 



2 fc afc t cd 



(11) 



III. 



SECOND ORDER SYMMETRY 
OPERATORS 



In this Section we shall solve (10). We begin by con- 
sidering the first equation therein. The coefficients E ab 
are zero-order matrix operators which can be expanded 
in the basis of C(2) defined above: 



E' 



ab „ab 



e ab I + ef 7 C + e ab 7 



(12) 



where the coefficients e ab , ef, e ab are point functions in 
M. Since the matrices I, j c , 7 are independent, the first 
equation in (10) implies 



e (ab e c)d = 
B (ab e c)d = q 



(13) 



By multiplying the second equation of (13) by e^e one 
obtains directly e^ ab 8 c ] = 0, the trace of which gives 
2e ab + e ab + e ab = 0. It follows immediately that 

e ab = (14) 

The first equation of (13) may now be solved to obtain 

ef 7 C = 2c*( a 7 b ) (15) 



Here a a are the coefficients along the frame e a ■= e^d^ 
of an arbitrary vector field a. 

Hence the most general E ab which satisfies the first 
equation in (10) is of the form 



E ab = e abj + 2a (« 7 fc) 



(16) 



We now consider the second equation of (10). By set- 
ting 



F a = / a I + / 6 a 7 b + / a 7 
it may be expanded to give 

'y( a a b ) = 

-2/ c (a e b ) c = -2V c a (a e b *> c 
_2/( a e b ) c = V d e ab 



(17) 



(18) 



The first of these equations are the Killing equations 
for the vector a a , which is hence a Killing vector or pos- 
sibly the zero vector. 

The second equation in (18) can be recast as 

(fc a — V c a( a ^ e b ' c = 0, from which one can easily prove 

that the most general 4>c satisfying the equation 4> c a e b "> c = 
is = AS®, where A is an arbitrary function on M. 
Hence we have 

f c = X7 c a a + A5 a c (19) 
The third equation in (10) can easily be rewritten as 
f a S b + f b 5 a c = e dc VV b (20) 
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Taking the trace we obtain 

f a = U bc V b e ac 



(21) 



which shows that f a is uniquely determined by e ac . How- 
ever, (20) contains six scalar equations, of which only one 
has been used. The other five equations are exploited by 
back-substituting (21) into (20) to obtain an equation for 
e ac alone, namely: 



e ef W e e af 5 b c + e ef V e e bf 8 a c = 3e dc V V b 



(22) 



This is an integrability condition for e ab which one can 
show is equivalent to 

V ( V b ' = (23) 

This equation implies that e ab is a Killing tensor which 
we shall from here on denote by K ab . We note that (23) 
is necessary condition for (18) to be true. 



To summarize the most general solution of the first and 
second equations in (10) is given by 



J E ab = K ab l + 2a( a 7 b) 

j F a = jaj + ( 7 c Vc a a + Aj a ) + \e bc V b K ac ~i 

where 



(24) 



( V (d K ab) = o 

jv< a a b ) = 



(25) 



We observe that there is yet no condition on the vector 
field f a . 

We now consider the third equation of (10). Writing 



G = gl + g a 7 a + .97 
we obtain the equations 



(26) 



V a A = =>ieC 

g b = -f a b (27) 
-2ge ab = V b f a + \e dc e e b V e V d K ac + §K ab + §K cd e d a e c b 



The first equation implies that A is constant, while the 
second uniquely determines g b in terms of a a and the 
Ricci scalar. To completely use the information con- 
tained in the third, which consists of nine scalar equa- 
tions, we consider separately its symmetric and antisym- 
metric parts. They may be written as 

= V( b / a ) + \e dc e e { N e \7 d K^ c + %K ab + §K cd e d {a e c ' 
4.9 - e ba W b .f a - \e dc V a V d K ac 

(28) 

Since K ac is a Killing tensor one can prove 
e dc V a 'V d K ac = 0. Hence the second equation in (28) 
implies that 

g = hbaV b r (29) 

Turning our attention to the first equation in (28), we 
rewrite it as 

V(« (/*>) - V c X b)c ) = (30) 

It follows that there exists a Killing vector £ a such that 
( a = f a — V c K ac from which we obtain 

f a = ( a + \7 c K ac (31) 

Hence the general solution of the first three equations 



in (10) is given by 
E ab = K ab l + 2a< a 7 b ) 

F a = + Vc x^)I + ( 7 c V c a a + Aj a ) + ±e bc V b K ac j 
G = gl- f a bl b + ie fca V"C a 7 

(32) 
where 

AeC 

V (d K ab) = o (33) 

V (a a b) = 0, V (d C h) = 

It remains to consider the fourth equation in (10). We 
would expect that this matrix equation would yield three 
equations, one each from I, 7 , and 7. However, in this 
case two of them are identically satisfied and the only 
additional condition is 

V a g = -\V b (RK a b ) ^d^g = -\V v (RK^) (34) 

This equation locally detemines g if (and only if) the 
right hand side is a closed 1-form. We may investigate 
the circumstances under which this condition holds with 
the use of some further information about Killing tensors 
on surfaces. For example, if M is flat the right hand side 
is trivially zero which implies that g € C is constant. 
However, if we consider the non-flat case, the condition 
is in general not satisfied. Its analysis requires knowledge 
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regarding which spin manifolds M admit non-trivial va- 
lence two Killing tensors. 

Fortunately, this knowledge is available. Indeed, one 
has the following result [3]: a 2- dimensional Rieman- 
nian space admits a non-trivial valence two Killing ten- 
sor if and only if it is a Liouville surface in which case 
there exists a system of coordinates (u, v) with respect to 
which the metric g and Killing tensor K have the follow- 
ing forms: 



g = (A(u) + B(v)) (du 2 + dv 2 ) 



(35) 



K = 



B(v) 



A(u) + B(v) 



d u ®d u - 



A(v) 



A(u) + B(v) 



d v ®d v = K ab e a ®e b 
(36) 



where A and B are arbitrary smooth functions. Further- 
more, the frame components of K are given by 



[K ab ] = diag{B{v), -A(u)) 



(37) 



Since in this case we know the form of the metric (and 
hence the scalar curvature R) and the Killing tensor, we 
may expand the integrability condition of equation (34) 
to obtain the following condition: 



(A + B) 2 {A'B"' + A"'B') + 6A'B' ((A') 2 + {B') 2 ) - 6A'B'{A + B)(A" + B") = 



(38) 



r 



where the primes denote differentiation with respect to 
the appropriate variable. 

If at least one of the two functions A and B is constant, 
then this integrability condition is trivially satisfied. In 
this case the Liouville surface M is a surface of revolu- 
tion, since it admits a Killing vector. For example if A 
is constant, then u is cyclic with d u the corresponding 
Killing vector. In this case, equation (34) determines g 
uniquely modulo a constant. 

To complete the analysis of condition (38) we need to 
investigate whether there exist solutions with neither^! 
nor B constant. Dividing by A'B' when needed and dif- 
ferentiating we obtain the following separated equation: 



(39) 



which we may easily integrate obtaining 

f (A') 2 = kA 4 + a 3 A 3 + a 2 A 2 + a x A + a 
I (B') 2 = -kB 4 + b 3 B 3 + b 2 B 2 + b x B + b Q 




(40) 



where k is the separatiuon constant, and ai and bi are 
integration constants. 

Substituting back for the second derivatives of (38) we 
obtain 

3(a 3 -b 3 )(A + B)=2(a 2 +b 2 ) (41) 
It follows that there are two cases to be considered: 

• Case I (a 3 ^ 63): then A and B must both be 
constant. This conclusion contradicts the original 
hypothesis that both A 1 and B' are non- vanishing. 



• Case II (ci3 — b 3 ): then a 2 = —b 2 . By substituting 
back into (38) we find that a\ = b\ and oq = — feo- 

We conclude that in the generic case at least one of 
the functions A and B must be constant. The only other 
possiblities are the special cases where A and B are so- 
lutions of the following system of differential equations: 



(A 1 ) 2 = kA 4 + a 3 A 3 + a 2 A 2 + a\A + a 
(B') 2 = -kB 4 + a 3 B 3 - a 2 B 2 + a x B - a Q 



(42) 



These special cases include, for example, the solution 



A{u) 



and B(v) = av which corresponds to the 



plane in parabolic coordinates. They also account for all 
constant curvature surfaces and some additional Liou- 
ville surfaces. In fact, one can compute the Ricci scalar 
for special cases (42) obtaining 



R= (A 



B)k + \a 3 



(43) 



We observe that by setting k = in (42) we obtain the 
surfaces of constant curvature which includes the plane 
when a 3 = 0. However, for k ^ we obtain more general 
surfaces. Which special class of Liouville surfaces satisfies 
condition (42) remains to be investigated. 

In the sequel, we shall consider the generic case in 
which either A' = or B' = 0. The details about special 
cases satisfying (42) will be considered in a forthcoming 
paper. 



IV. FIRST ORDER OPERATORS AND 
REDUCTIBILITY OF SECOND ORDER 
OPERATORS 



First order operators are obtained by setting 



and a a = 0. The integrability condition is trivially 
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satisfied and g G C. Thus the most general first order 
operator has the form 



'E afc = 

G = gl+\e ba V b C a 1 



where 



A, g e C 



(44) 

Among first order operators we identify the trivial 
ones, namely the first order operators of the form 



aB + /3K 



(45) 



for some zero order operator Ko = gl. The following 
theorem characterizes the trivial first order operators: 

THEOREM: There is a one-to-one correspondence be- 
tween non-trivial first-order symmetry operators and 
Killing vectors £ on M. ■ 

The situation for second order symmetry operators is 
somewhat more complicated although it can be analyzed 
along similar lines. 

We define the trivial second order symmetry operators 
as the ones of the form D o IK^ + = K4 o D + K.[ , where 
Ki and denote any first order symmetry operators. 
Thus the most general non-trivial second order symmetry 
operator has the form 



where 



E ab = K ab l 

¥ a = V c K ac l + \e bc V b K ac -y 
G = gl 

R ab Xr] ab 

V[, (V* (RK v]x ))=0 



(46) 



(47) 



and where g a solution of equation (34) . This results is 
summarized in the following theorem: 

THEOREM: There is a one-to-one correspondence be- 
tween non-trivial second order symmetry operators and 
Killing tensors K ab on M such that K ab ^ Xr] ab and 
V [fl {V x (RK„ ]x )) =0. 

In particular, if we have a non-trivial second order sym- 
metry operator on M, then M must admit a Killing ten- 
sor K ab ; hence M is a Liouville surface (see [3]). Ac- 
cording to the analysis done in Section 3, M generically 
must admit a Killing vector £ (corresponding to the cases 
A' = or B' = 0) or, non-generically, it must be in the 
form (42). 

A reducible second order symmetry operator is a linear 
combination of products of first order operators, namely 
A* J o Kj 1 for any first order symmetry operators Kj . 

If we consider the product of two first order operators 
with = C(")I + A (i)l a (i = 1) 2 ), w e obtain the second 
order operator defined by 



K a 



Si S2 



2A 1 A 2V ab a a = t(A 1 Q + A 2 (48) 



One can easily check that, since ^ are Killing vectors, 
K ab is automatically a Killing tensor and that a a is a 
Killing vector, since it is a linear combination of Killing 
vectors. 

On particular surfaces (for example, surfaces of con- 
stant curvature) all Killing tensors are known to be pro- 
duced by linear combinations of symmetrized products 
of Killing vectors. Thus on such surfaces all non-trivial 
second order symmetry operators are reducible. To have 
non-trivial irreducible second order symmetry operators 
one needs Killing tensors other than the ones which are 
symmetrized products of Killing vectors. 

The conditions for the existence of irreducible Killing 
tensors can be shown to be met on a explicit example. 
Let M e an asymmetric ellipsoid endowed with the metric 
induced by that of the Euclidean Space. As shown by Ja- 
cobi, the Hamilton- Jacobi equation of the geodesies of M e 
admits orthogonal separable coordinates. Therefore, also 
a symmetric Killing 2-tensor K e associated with them ex- 
ists (see for example [1]). However, no Killing vector can 
exist on M e . It follows that K e is irreducible. 

Indeed, let us consider 



-b 2 



1 



(49) 



with < a 2 < u\ < b 2 < u\ < c 2 < u% and (x, y, z) are 
Cartesian coordinates. The equations (49) for i = 1 ... 3 
define three families of confocal quadrics parametrized by 
Ui] u\ and u 2 associated with one and two- folded hypcr- 
boloids, u 3 with the ellipsoids. By defining M e as u 3 = h 
for some constant h, the variables U\ and u 2 determine 
orthogonal separable coordinates on M e . From (49) it is 
possible to write 



_ {a 2 -h 2 )(a 2 -ul)(a 2 ~ul) 



X" = — 
„,2 _ (b 2 



(a 2 -b 2 )(a 2 - 
-fe 2 )(b 2 -^)(b 2 



-ul) 



(a 2 -b 2 )(b 2 -c 2 ) 
(c 2 -h^-ul){c 2 -ul) 
(a 2 -c 2 )(b 2 -c 2 ) 



(50) 



and to obtain 



■ yn ~ (u 2 -c 2 )(b 2 -u{)(a 2 -u 2 ) 
u\(u\-h 2 )(u\-u\) 
(u 2 -c 2 )(b 2 -u 2 )(a 2 -u 2 ) 



922 



(51) 



It is evident that a rescaling u = u(ui), v — v(u 2 ) of the 
coordinates puts the metric in Liouville form 



g vv = A(u) + B(v) 



(52) 



with the two functions A and B both non-constant. 
Then, the coordinates (u, v) are separable and the associ- 
ated Killing tensor is given by (36). Such a Killing tensor 
is of course irreducible since the metric (51) admits no 
Killing vector. We remark that the Killing tensor given 
in this case by (36) does not satisfy condition (42); This 
is consistent with the fact that Dirac equation does not 
separate in this coordinate web (see [15]). Whether this 
feature, namely that existence of Killing vectors, is a nec- 
essary condition for the separation of the Dirac equation 
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and for the existence of second order symmetry operators 
also in the special cases (42) in addition to the generic 
case, will be investigated in a forthcoming paper. 



V. SECOND ORDER SYMMETRY OPERATORS 
AND SEPARATION OF VARIABLES 

It is known ([20], [11] and references therein) that 
the Dirac equation separates only in cooordinate sys- 
tems separating the Hamilton- Jacobi equation for the 
corresponding metric tensor, and only if one at least of 
these coordinates is associated to a Killing vector. In 
2-dimensions, this is the case for the coordinate system 
(it, v) for which the metric is in Liouville form, as seen in 
Section 3, where, say, u is an ignorable coordinate. The 
appropriate spin-frame in this case is given by 







{<) = 



y/A(u)+B(v) 





(53) 



where A = and B = (3(v)~ 2 . Assuming that vectors 
a and £ are zero the Dirac operator B and second order 
symmetry operator K may be written in matrix form as 



-1 

1 



d u 



i 
-i 



8y 



-13' 
/3' 



where primes denote differentiation and 



(54) 



(55) 



These expressions for E) and K coincide with the separa- 
tion scheme D 5 and the corresponding symmetry oper- 
ator L 5 found in [15] 1 obtained by analyzing naive sep- 
aration for the Dirac equation in two-dimensional Rie- 
mannian manifolds. It is an instance of non-factorizablc 
separation [12]. By computing 



mip = 



and 



Kip - fiijj = 
with the separation ansatz 



ai(u)bi(v) 
a 2 (u)b 2 (v) 



(56) 



(57) 



(58) 



1 In [15], the right-hand side of each equation expressing quantities 
X and Y as functions of A;, B{ etc., from Section 7 on, must 
be divided by 2; this correction is necessary to compare with 
notations used here and it affects all expressions derived from 
the formulae of above. In particular, due to additional typos, the 
right-hand term of the last equation in (7.3) should be multiplied 
by i k/2. 



aft 



;er dividing the two components of Hip by j3a\b 2 and 
f3a 2 b\ respectively, we obtain the separated equations 



..&' 



'hi 



" = j 12. — iSL hs. 4- m in 

.02 bi 2/3 61 "T 61 



2/3 b 2 



m 62 



m 61 
b 2 



(59) 



and, from (57) we have 




(60) 



By introducing separation constants \i\ and [i 2 such that 
M1M2 = Mi the first system gives for the aj 



and for the h; 



a 2 = -fiitti 
a[ = \x 2 a 2 



■28 



ib' 2 - i§pb 2 + r -%b 2 = n 2 b\ 



(61) 



(62) 



It is evident that K provides, via (60), a decoupling re- 
lation for the equations (61). On the other hand, (60) 
can be obtained by applying twice equations (61). The 
symmetry operators associated with separation are gen- 
erated by this means in ([15]). The solutions for the a, 
are 




ci sin(y^*u) + c 2 cos(^fjiu) 

c 2 sin(^/7iu) + a cos(^/Jlu)] 



(63) 



where c\,c 2 € C. 

In general, the relations (62), cannot be associated 
with symmetry operators (they can, for example, if /3 
is a constant, that is for Cartesian coordinates). How- 
ever, again it is possible to use them to obtain decoupled 
second-order differential equations in the 6j and then the 
solutions of (62) itself, in the same way as for the dj. The 
solution obtained in this way is consistent; indeed, the 
separation constant appearing here is again /i = /iifx 2 
and it is the eigenvalue of the symmetry operator IK. 
Therefore, no additional symmetry operator is necessary 
for the separation. The general solution of (62) can be 
easily computed in Cartesian coordinates, (3(v) — 1, and 
we have 

\b 1 =d 1 sin(Mw) + d 2 cos(Mv) 

\b 2 = -^[{di +iMd 2 )sm(Mv) + (d 2 - iMdi) cos(Mv)] 

(64) 

where d\ , d 2 G C and M = \Jm? — \i. It is remarkable 
that, even if the at and the 6j depend on the fij, the 
products ipi = ciibi depend on /i only. 

By setting, for example, f3 equal to 1, e v , sinh(w), 
cosh(w), k — cos(v), with k > 1, respectively, the cor- 
responding Riemannian manifolds are: the Euclidean 
plane in Cartesian and polar coordinates, the sphere, the 
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pseudo-sphere (or hyperbolic plane), the torus with met- 
ric induced from that of the Euclidean Space. 

For polar coordinates in the Euclidean plane and in 
the cases of sphere and pseudosphcrc, equations (62) can 
be (formally) integrated with the use of computer alge- 
bra (Maple), to obtain solutions respectively in terms of 
Bessel functions and, on the sphere and pseudo-sphere, 
in terms of hypergeometric functions. 

To summarize, we obtained various examples of differ- 
ent behaviours with respect to Dirac equation separabil- 
ity. The skew ellipsoid is an example in which one has an 
(irreducible) Killing tensor and separable coordinates for 
Hamilton- Jacobi equation though this Killing tensor is 
not associated to a second order symmetry operator for 
the Dirac equation; this is consistent with the fact that 
Dirac does not separate in these coordinates. In the case 
of parabolic coordinates one has a Killing tensor, sep- 
arable coordinates for Hamilton- Jacobi equation and a 
second order symmetry operator for Dirac, although the 
Dirac equation does not separable in these coordinates. 
Finally, we presented several examples for which one has 
a Killing tensor, separable coordinates for the Hamilton- 
Jacobi equation and second order symmetry operators 
for which the Dirac equation is found to separate in the 
same coordinates. 



VI. CONCLUSION 

We compute the most general second order symme- 
try operator for the Dirac equation on a 2-dimensional 
spin manilfold. Particular choices of coodinates, spin 



frames or Dirac matrices play no role in the calculations 
which is manifestly covariant. The case of 2-dimensional 
Lorentzian spin manifolds will be considered in a subse- 
quent paper. We hope to use the insight obtained in this 
paper and in [15] to give an invariant characterization 
of separability for Dirac equation comparable to that for 
the Schrodinger equation described, for example, in [13]. 
It is unclear whether it will be possible to extend the 
same approach to higher dimensions, although we hope 
to be able to find new theoretical tools to tackle the prob- 
lem. The expectation is to be able to study the orbit of 
the group of spin transformations, which acts canonically 
on symmetry operators and consequently on the under- 
lying separation schemes. This information is necessary 
in view of a classification of non- equivalent separation 
schemes. In fact, the action of spin transformation group 
defines an equivalence relation among symmetry opera- 
tors which may define non-trivial equivalence between 
different schemes. 
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